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We study the thermodynamics of the relativistic quantum O(TV) model in two space dimensions. 
In the vicinity of the zero-temperature quantum critical point (QCP), the pressure can be written 
in the scaling form P(T) = P(0) + A r (T 3 /c 2 )J r j V (A/T) where c is the velocity of the excitations at 
the QCP and A is a characteristic zero-temperature energy scale. Using both a large-TV approach 
to leading order and the nonperturbative renormalization group, we compute the universal scaling 
function J-jv. For small values of TV (TV < 10) we find that Tn^x) is nonmonotonous in the quantum 
critical regime < 1) with a maximum near x — 0. The large- TV approach - if properly interpreted 
- is a good approximation both in the renormalized classical (x < —1) and quantum disordered 
(x > 1) regimes, but fails to describe the nonmonotonous behavior of J-jv in the quantum critical 
regime. We discuss the renormalization-group flows in the various regimes near the QCP and make 
the connection with the quantum nonlinear sigma model in the renormalized classical regime. We 
compute the Berezinskii-Kosterlitz-Thouless transition temperature in the quantum 0(2) model and 
find that in the vicinity of the QCP the universal ratio Tbkt/Ps(0) is very close to it/2, implying that 
the stiffness p„(T^ KT ) at the transition is only slightly reduced with respect to the zero-temperature 
stiffness p s (0). Finally, we briefly discuss the experimental determination of the universal function 
Ti from the pressure of a Bose gas in an optical lattice near the superfluid-Mott-insulator transition. 
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I. INTRODUCTION 

Many zero-temperature critical points observed in 
quantum many-body systems are described by a relativis- 
tic effective field theory [1, 2]. Bosonic cold atomic gases 
constitute a very clean experimental realization of such 
quantum critical points (QCP): a Bose gas in an optical 
lattice undergoes a quantum phase transition between a 
Mott insulator and a superfluid state [3-6]. When the 
transition occurs at fixed density, it is described by a 
relativistic quantum 0(2) model [7, 8]. 

Recent works have focused on the excitation spectrum 
of the relativistic quantum O(A) model in the vicinity 
of the QCP and in particular on the spectral function of 
the amplitude ("Higgs") mode in the broken-symmetry 
phase [2, 9-12]. Signatures of the amplitude mode have 
recently been observed in a two-dimensional superfluid 
near the superfluid-Mott-insulator transition [13]. 

In this paper, we study the thermodynamics of the rel- 
ativistic quantum O(A) model in two space dimensions. 
We extend previous results [1, 14] obtained to leading 
order in the large-TV limit by computing the full scal- 
ing function J-^ (x) determining the temperature depen- 
dence of the pressure near the QCP. Using a nonpertur- 
bative renormalization-group (NPRG) approach [15-17], 
we then calculate ^^(x) for finite values of A, including 
TV = 2 and A = 3. 



We start from the action 

where we use the shorthand notation 

X = (r,r), J dx = J dr J d 2 r. (2) 

<p{x) is an TV-component real field and r g [0, 0\ an imag- 
inary time (j3 = 1/T and we set h = ks = 1). and uq 
are temperature-independent coupling constants and cq 
is the (bare) velocity of the (p field. The factor 1/TV 
in Eq. (1) is introduced to obtain a meaningful limit 
TV — > oo (with uq fixed). The model is regularized by 
a ultraviolet cutoff A. In order to maintain the Lorentz 
invariance of the action (1) at zero temperature, it is 
natural to implement a cutoff on both momenta and fre- 
quencies but we will also sometimes use a cutoff acting 
only on momenta. 

In two space dimensions, the phase diagram of the rel- 
ativistic quantum 0(A) model is well known (Fig. 1) [ ]. 
At zero temperature, there is a quantum phase transi- 
tion between a disordered phase (ro > tq c ) and an or- 
dered phase (ro < r$ c ) where the 0(A) symmetry of the 
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FIG. 1. Phase diagram of the relativistic 0(N) model in two 
space dimensions for N > 3 [Eq. (1)]. The thick line shows the 
zero-temperature ordered phase with long-range order (LRO), 
while the dashed lines are crossover lines between the renor- 
malized classical (RC), quantum critical (QC) and quantum 
disordered (QD) regimes. The dotted line shows the limit 
of the high-T region where the physics is not controlled by 
the QCP anymore. (For N — 2, there is a finite-temperature 
BKT transition line for ro < ro c , which terminates at T — 
for r = r 0c .) 



action (1) is spontaneously broken (uq and Co are con- 
sidered as fixed parameters). The QCP at r = r 0c is in 
the universality class of the three-dimensional classical 
0(A) model with a dynamical critical exponent z = 1; 
the phase transition is governed by the three-dimensional 
Wilson-Fisher fixed point. At finite temperatures, the 
system is always disordered for N > 2, in agreement 
with the Mermin- Wagner theorem, but it is possible to 
distinguish three regimes in the vicinity of the QCP: a 
renormalized classical regime, a quantum critical regime, 
and a quantum disordered regime [1, 18]. For N = 2 
and ro < ro c , there is a finite-temperature Berezinskii- 
Kosterlitz-Thouless (BKT) phase transition [19-21] and 
the system exhibits algebraic order at low temperatures. 
The BKT transition temperature line Tbkt terminates 
at the QCP r — r 0c . 

Below the upper critical dimension d+ = 3 (d+ + z = 4) 
of the quantum phase transition, we expect the hyper- 
scaling hypothesis to hold. In two dimensions, this allows 
us to write the pressure in the critical regime as [22] 



P(T) = P(0) + N%F N 



T 



(3) 



where Fn is a universal scaling function, c the velocity 
of the critical fluctuations at the QCP [23] and |A| = 
| A(ro)| a characteristic energy scale at zero temperature. 
When r > r 0c , the system is disordered and we choose A 
to be equal to the excitation gap m$ oc (r — r 0c ) ziy of the 
ip field [y denotes the correlation-length exponent at the 
QCP) - not to be confused with the amplitude ("Higgs") 
mode gap. When ro < tq c it is convenient to take A neg- 
ative such that — A is the excitation gap in the disordered 
phase at the point located symmetrically with respect to 
the QCP, i.e. |A(r )| = m (2r 0c - r ) [ ]. -A is then 
proportional to the stiffness p s , the ratio \A\/p s being 
universal. With these definitions, A varies from nega- 
tive to positive values as we go across the QCP coming 



from the ordered phase. The two crossover lines shown 
in Fig. 1 are roughly defined by |A| ~ T. We stress that 
the scaling function F N is independent of all microscopic 
parameters of the model such as ro, uq or cq. The latter 
enter the temperature variation of the pressure [Eq. (3)] 
only indirectly via the values of the renormalized velocity 
c and the energy scale A. 

In the critical regime near the QCP, all thermodynamic 
quantities can be written in a scaling form. In addition 
to Fn, we will compute the universal scaling function Pjv 
which determines the excitation gap 



m(T) = TP 



N 



(4) 



at finite temperatures. As we shall see, the knowledge of 
Fn is necessary to obtain Fn in the large- TV limit. 

The outline of the paper is as follows. In Sec. II, we 
compute the universal scaling functions Fn and Fn to 
leading order in a 1/N expansion. We then use a NPRG 
approach to calculate F N and Fn for any value N > 2 
(Sec. III). The main results are presented in Sec. IIIB. 
Section III C is devoted to a detailed analysis of the RC 
flows in the renormalized classical, quantum disordered 
and quantum critical regimes for N > 3. In the renor- 
malized classical regime, where the physics is dominated 
by the N — 1 Goldstone modes of the zero-temperature 
broken-symmetry phase, we show that the NPRG flow 
equations yield the one-loop RG equations of the quan- 
tum O(N) nonlinear a model (NLerM) [ ]. The BKT 
transition temperature in the quantum 0(2) model is dis- 
cussed in Sec. Ill E. The implication of our results for cold 
atomic gases are briefly discussed in the Conclusion. 



II. LARGE- iV LIMIT 

In this section, we use a cutoff A acting only on mo- 
menta, i.e. |q| < A. We do not distinguish between the 
bare velocity cq and the renormalized one c since they 
coincide in the large- A limit. 

Following the standard method in the large-A limit 
(see, e.g., Refs. [24, 25]), we express the partition function 
as 



Z = J T>[(p, p, A] cxp|— J dx 
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(5) 



It can be easily verified that by integrating out A and 
then p, one recovers the original action S[(p\. If, instead, 
we first integrate out p, we obtain 



V[ip, A] exp 



dx 



2 (V ^ + 2>^ 2 



• A 2 
■*2* 



3JV / , /-A N2 

+ — /^A-ro) 2 



(6) 
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We then split the ip field into a field a and a (N — 1)- 
component field 7r. The integration over the 7r field gives 

j p[ 7r ] e -/ d -[i( v -) 2 +^(^-) 2 +^- 2 ] = (det.9)^- 1 )/ 2 , 

(7) 

where 

g~ x {x, x') = [- V 2 - c7 2 d 2 + i\(x)]5(x - x') (8) 

is the inverse propagator of the iri field in the fluctuating 
A field. We thus obtain the action 

S[a, X] = ^Jdx [(V<t) 2 + c~ 2 (d r o) 2 + iXa 2 ] 



3N 
2u Q 



dx (iX — ro) 2 + 



N — 1 



Trln g 



(9) 



In the limit N — > oo, the action becomes proportional 
to N (this is easily seen by rescaling the a field, a — > 
VN a) and the saddle-point approximation becomes ex- 
act. For uniform and time-independent fields <j{x) = a 
and X(x) = A, the saddle-point action is given by 

w SM= >'-^- T ' r+ m^ (10) 

(we use iV — 1 ~ N for large N), with g^ 1 (q) = q 2 + 
w 2 /c 2 + iX in Fourier space, q — (q, ioj n ), u> n = 2-kTu (n 
integer) is a bosonic Matsubara frequency, and V denotes 
the volume of the system. From (10), we deduce the 
saddle-point equations 



am 1 = 

9 6A / m 2 

a = ~ 

u \ c 1 

where we use the notation 



ro 



N I g(q), 
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(11) 



d 2 q 
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where f = f°° ^. 

In the disordered phase ro > r$ c , a = and the mass 
rn = m(T = 0) is determined by Eqs. (11,13), 
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j 2 + c 2 q 2 + TOq uj 2 + c 2 q 2 



6toq too 6 

M C 2 47T U 



{ro - roc). 



(14) 
= 0, (15) 

(16) 



By comparing the two terms in the lhs of this equation, 
we obtain a characteristic momentum scale, the Ginzburg 
scale ka ~ cuq/24it, which signals the onset of critical 
fluctuations [ ]. In the critical regime, mo <C cfcc, we 
obtain 



2A7T 

m = (ro - r 0c ), 

u 



(17) 



which gives zv = 1, i.e. a correlation-length exponent 
v = 1 since the dynamical critical exponent z = 1. In 
the noncritical regime too 3> cfcfj, TO o ~ ( r o — ^Oc) 1 ^ 2 
and we recover the classical value v = 1/2. The anoma- 
lous dimension rj vanishes to leading order in the large- N 
limit. 

In the ordered phase ro < ro c , mo vanishes and a is 
finite, 



-6N— - N 



U 



q J u 



lj 2 + c?q^ 



(ro -r 0c ). 



(18) 



(12) The stiffness is equal to p s = Na 2 [27]. 



and m 2 = iAc 2 (iX is real at the saddle point). These 
equations show that the component a of the ip field which 
was singled out plays the role of an order parameter. In 
the ordered phase, a is nonzero and m = 0. The propa- 
gator g(q) — l/(q 2 +cj 2 /c 2 ) is gapless, thus identifying 
the TTi fields as the N — 1 Goldstone modes associated 
with the spontaneously broken O(N) symmetry. In the 
disordered phase, a vanishes and m determines the gap 
(or "mass" ) of the cp field as well as the correlation length 
£ = c/m. 



A. Zero temperature 

The critical value ro c corresponding to the QCP sepa- 
rating the ordered and disordered phases is obtained by 
setting fj — to == in Eqs. (11), 
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(13) 



B. Finite temperatures 

At finite temperatures, the system is always disordered 
((7 = 0), in agreement with the Mermin- Wagner theorem, 
and the mass m is obtained from the saddle-point equa- 
tion 



o = A 

6 

u 
G 



m 

^2 



m 



ro 



'•() 



c 2 q 2 



T ( sinh ^ 
2^ n l^sinhf, 



Uo 



i \ T ( ^ ■ m \ 6to 2 

(roc - r ) + — In ( 2 sinh — ) + 



2tt 



2TJ 



(19) 



In the critical regime the last term can be neglected and 
we obtain 



= 2T asinh 



(20) 
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where we have introduced the characteristic energy scale 
|A| defined by 



a 247r , X 

A = (r - r 0c ). 

u 



(21) 



A corresponds to the T = gap m on the disordered 
side r > r 0c of the QCP, and to —4irp s /N on the ordered 
side ro < r$ c with p s the zero-temperature stiffness (see 
the discussion in the Introduction). The critical regime 
is defined by T, |A| < ck G . 
We can rewrite Eq. (20) as 



with the universal scaling function 



Foo{x) = 2asinh -e 



1 



r/2 



(22) 



(23) 



-Foo^) satisfies 



Foo(x) 



x/2 



if a; 



2asinh(l/2) if x 
x if x - 



—00, 
: 0, 
00, 



(24) 



with 2asinh(l/2) ~ 0.962424. The three cases in 
Eq. (24) correspond to the renormalized classical (m ~ 
Te -\A\/2T^ quantum cr itical (m ~ 2 asinh(l/2)T), and 
quantum disordered (m ~ A) regimes, respectively (see 
Fig. 1). 



C. Pressure 

In the large- TV limit, the pressure P = —S[a, X]//3V is 
obtained from the saddle-point value of the action, 



P 

N 



2Nc 2 2u Q 



ro 



1 



2f3V 



Trine/ -1 . (25) 



Using the results of Appendix A for Trlng -1 , in the 
critical regime we can write the pressure in the scaling 
form (3) with the universal scaling function 
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Foo{x)lA2{t 



Li s (z) is a polylogarithm, 



fe=i 



(26) 



(27) 



and O(x) denotes the step function. 



From the definition of Foo{x), we obtain the limiting 
cases 



•F 30 \X) 



03) 
2tt, 



if x — > —00, 



2C(3) s 



if x = 0, 

t)7T 

if x — > 00, 



(28) 



where Ci z ) is the Riemann zeta function: £(3)/27r ~ 
0.191313 and 2C(3)/5tt ~ 0.153051. To obtain 
linLj._j._oo Foo(x) we use F 00 (x) — > e x l 2 for x — > —00 and 
I_3(l) = C(3). The universal number Foo(0) is obtained 
noting that (0) = 2 In r, with r = (1+-/5) /2 = 2-t~ 2 
the Golden mean, and using [28] 



1 



Li 2 (2-r) = — - -ln 2 (2-r) 
; 15 4 v y 



Li 3 (2-r) = -C(3) 



^M2-r)-lln 3 ( 2-r). 



(29) 



It should be noted that the scaling function Foo(x) 
as well as -^(O) agree with results obtained from the 
NLcrM in the large-- / limit [1, 14, 28]. This follows from 
the fact that the linear and nonlinear O(A) models are 
in the same universality class and therefore exhibit the 
same critical physics. 



III. NPRG APPROACH 

The strategy of the NPRG is to build a family of the- 
ories indexed by a momentum scale k such that fluctu- 
ations are smoothly taken into account as k is lowered 
from the microscopic scale A down to [15-17]. This is 
achieved by adding to the action (1) the infrared regula- 
tor 



&S k [<p} = \ ^^i{-q)Rk{q)fi{q) 1 



(30) 



so that the partition function 

Z k [J] = J T>[<p]e- s M- ASk M + f dx ^ JiVi (31) 
becomes k dependent. The fc-dependent effective action 
T fc [4>] = - In Z k [J] + I dx Ji<t>i - [0] (32) 

is defined as a modified Legendre transform of — In Zf~ [J] 
which includes the subtraction of ASk[<p]- Here <p(x) = 
(ip(x)) is the order parameter (in presence of the exter- 
nal source). The initial condition of the flow is speci- 
fied by the microscopic scale k = A where we assume 
that the fluctuations are completely frozen by the ASk 
term, so that I\[0] = S[<f>]. The effective action of the 
original model (1) is given by Tk=o provided that Rk=o 
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vanishes. For a generic value of k, the cutoff function 
Rk(q) suppresses fluctuations with momentum |q| < k 
or frequency \uj n \ < c k k but leaves unaffected those with 
|q|,|wn|/cfc > k (here Cfc denotes the (renormalized) ve- 
locity of the cp field). The variation of the effective action 
with k is given by Wetterich's equation [ ] 



d t r k [<p] = ~ti^ v?/, 



(2) 

where t = \n(k/A). T k [<f>] denotes the second-order 
functional derivative of Tj. [</>]. In Fourier space, the trace 
involves a sum over momenta and Matsubara frequencies 
as well as the internal index of the <f> field. We use a reg- 
ulator function R k {q) which acts both on momenta and 
frequencies, 



(33) 



Rk(q) = Z A ,k q 



,2 \ / „2 

r 



.,:-\ itf+ui/ci 



k 2 



(34) 



where r(Y) = l/(e Y — 1). The fc-dependent constant 
Z Aik is defined below [Eq. (37)]. 

When <j> is constant, i.e. uniform and time indepen- 
dent, the effective action coincides with the effective po- 
tential. 



U k {p) 



1 

W 



T k [4>] 



<p const 



(35) 



Because of the O(N) symmetry of the effective action Tk, 
the effective potential Uk{p) must be a function of the 
O(iV) invariant p = <f> 2 /2. The pressure is then simply 
defined by 



P(T) = -U k = ( Po ), 

where po,fe denotes the position of the minimum 
and p = limfc-M) p ,k- 



(36) 
of U k (p) 



A. Approximate solution of the flow equation 



Because of the regulator term ASk, the vertices 

k,ii ■••i n 



LPA, and therefore the LPA', are exact in the large-A^ 
limit [30]. To further simplify the analysis, we expand 
Uk{p) about the position p 0) fc of its minimum, 



_ / U k (p Q>k ) + ^(p- Po,k) 2 if Po,k > 0, 



Uk{p) -^ UkM + 5kP+^p 2 



if Po,, 



0. 



(38) 



Although the RG equations can also be solved for the 
full effective potential, the determination of the singular 
part of the pressure turns out to be extremely difficult in 
that case due to numerical noise. 

The LPA is known to be very accurate to obtain ther- 
modynamic quantities. It has been used to compute 
the pressure in the three-dimensional quantum ip 4 the- 
ory with Ising symmetry (i.e. N = 1) [31, 32]. The 
results compare very well with those of the Blaizot- 
Mendez-Wschebor approach (BMW) - an elaborated 
NPRG scheme which preserves the full momentum and 
frequency dependence of the propagator [33-35]. There 
are also strong indications that the LPA (or the LPA') is 
a good approximation even when it is supplemented by 
a truncation of the effective potential [Eq. (38)]. As will 
be shown below, the truncated LPA' remains accurate 

- and nearly exact in the renormalized classical regime 

- in the limit N — > oo [ ]. Furthermore, it has also 
been used to determine the phase diagram of the Bose- 
Hubbard model in two and three dimensions [8, 37, 38]: 
although a truncation of the effective potential leads to a 
loss of accuracy, the results remain within 10 percent of 
the exact ones obtained by quantum Monte Carlo simu- 
lation [39, 40]. 

The derivation of the flow equation for U k (p), Z Atk 
and V A ,k is standard [15, 16] (the only difference with 
the classical 0{N) model comes from the finite size f3 
in the imaginary-time direction [41, 42]). The effective 
potential satisfies the flow equation 

dtUk(p) = \J d t R k (q)[G k A(q;p) + (N-l)G k Aq;p)l 

(39) 

where 

G-]( g; p) = Z A ^ 2 + V A W n + U' k { P ) + 2pU k '(p) + R k (q), 



(<7i, • • • , q n ) are smooth functions of momenta G k \(q; p) = Z A ^ k q 2 + V A .kU 2 n + U' k {p) + R k {q) 



and frequencies and can be expanded in powers of qf/fc 2 
and w 2 . /c 2 k 2 . Thus if we are interested only in the long- 
distance (critical) physics, we can use a derivative expan- 
sion of the effective action [15, 16]. In the following, we 
consider the ansatz 



(40) 



determine the longitudinal and transverse parts of the 
propagator Gk — (T^ + Rk)^ 1 in a constant (i.e. uni- 
form and time-independent) field </>, 



Z A . 



■(V0) 2 



v A „ 



id T <p) 2 + u k {p) 



(37) 

which is often referred to as the LPA'. It differs from 
the local potential approximation (LPA) by the intro- 
duction of two field renormalization constants Z A _ k and 

V A}k (Z a ,a = 1 and Va,a = c q 2 )- It is the minimal 
ansatz beyond the LPA which includes a finite anoma- 
lous dimension rj at the QCP (see below). Moreover the 



Gk,ij(q;4>) = ~^~ Gk ^P) + ( S *,3 ~ 



2p 



G k , t (q;p). 
(41) 

The contribution of G k t to d t Uk comes with a factor 
N — 1 corresponding to the number of transverse modes. 
When po.fc > 0, U' k (po,k) vanishes and these modes be- 
come gapless for R k (q) — > (Goldstone modes). The 
stiffness is given by p Sj fc = 2Z A k po,k [27]. In the disor- 
dered phase, the minimum of U k (p) is located at po. k = 



so that all modes exhibit a gap m k = \/U' k (0)/VA yk (for 
Rk(q) 0) corresponding to a finite correlation length 
S,k — c k I m k where c k = U Za ^/Va^ is the renormalized 
velocity (ca = Co). The actual gap m and correlation 
length £ in the disordered phase are obtained for k = 0. 

The flow equations for ZA,k and Va,u are obtained from 
the flow equation (33) by noting that 
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Z AM = lim^rg(g;p ,fe), 



(42) 



At the zero-temperature QCP, ~ fc and V^fc ~ 

fc— i7-2(«-x) [ 4; j j ^ich allows us to deduce both the 
anomalous dimension r\ and the dynamical critical ex- 
ponent z. The latter is equal to one due to the Lorentz 
invariance of the action (1) at T = 0. The exponent v can 
be obtained from the divergence of the correlation length 
£ ~ ( r o — r Qc)~ v in the disordered phase as the QCP is 
approached, or more directly from the escape rate from 
the fixed point when the system is nearly critical. 
The RG equations are given by [44] 



d t po,k 
d t S k 



- » h,t if Po,fe > 0, 



At 



(N + 2)4,! if p o ,fe = 0. 



d t Xk = -A^[9J M (0) + (N - l)J fe)tt (0)], 


\Z A ,k = -2\ 2 k p 0:k — {J kM ( P ) + Jk,n(p)]\ p=0 , 



(43) 



dp 

d 



d t V A ,k = -2A fc p . fe7 — [Jfc.ti(p) + Jfc.it (p) 
out 



lp=0' 



while the equation for the thermodynamic potential (per 
unit volume) U k (po,k) is directly obtained from (39). We 
have introduced the threshold functions 

Ik,a= j d t G kt a{q\ po,k), 

r . ( 44 ) 
Jk, a p{p) = / [d t G k ^ a (q; p 0yk )]G k ^(p + q; po, k ), 

Jq 

with a,P = 1, t. The operator 3 t = (dtR k )da k acts only 
on the t dependence of the cutoff function R k . The prop- 
agators G k j(p;p Q . k ) and G k<t (p; p , k ) are given by (40) 
with U' k (pa,k) = fa and U'^(p ,k) = Afc. 

The flow equations are solved numerically. Results re- 
lated to the thermodynamics are discussed in the follow- 
ing section. 



B. Universal scaling functions 
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FIG. 2. Universal scaling function Fz(x) [Eq. (4)] computed 
for various values of the microscopic parameters (A = 100 and 
co = 1). 
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FIG. 3. Same as Fig. 2 but for the universal scaling function 
T 3 (x) [Eq. (3)]. 
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FIG. 4. Universal scaling functions Fn and J-jv for various 
values of N obtained from the NPRG. The black points show 
the analytic results (23) and (26) in the limit N — t oo. 



We first solve the equations at T — to determine 
roc and obtain the critical exponents v and r\ as well as 
the characteristic energy scale A = A(ro). For N = 3 
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N 1000 10 8 6 4 3 2 
p a /(JV|A|) 0.0838 0.0853 0.0864 0.0891 0.0965 0.1059 0.1321 



TABLE I. Universal ratio p s /(N\A\) in the zero-temperature 
ordered phase. The exact result in the limit N — > oo is 1/4tt ~ 
0.080. 



we find v ~ 0.699 and r\ = 0.0507, to be compared 
with the best estimates for the three-dimensional 0(3) 
model obtained from resummed perturbative calcula- 
tions [45] (y ~ 0.7060, rj = 0.0333), Monte Carlo sim- 
ulations [46] {v ~ 0.7112, r/ = 0.0375), or the NPRG in 
the BMW approximation [35] iy ~ 0.715, T] = 0.040). 
For N = 2, our results v ~ 0.613 and r? = 0.0582 
should be compared with the critical exponents of the 
three-dimensional 0(2) model: resummed perturbative 
calculations [ ] [y ~ 0.6700, rj = 0.0334), Monte Carlo 
simulations [47] (y ~ 0.6717, r\ = 0.0381), NPRG- 
BMW [35] (v ~ 0.674, 77 = 0.041). Note that the 
rather poor estimate of rj is a well-known limitation of 
the LPA'; a much better result can be obtained by con- 
sidering the full derivation expansion to order 0(d 2 ) [48]. 
At finite temperatures, the two-dimensional relativistic 
0(2) model exhibits a BKT phase transition. Although, 
stricto sensu, the NPRG does not capture this transi- 
tion, most universal properties of the latter are neverthe- 
less correctly reproduced [49, 50]. In particular, recent 
work on the two-dimensional Bose gas has shown that 
the thermodynamics can be accurately computed using 
the NPRG [51]. The BKT transition is further discussed 
in Sec. HIE. 

Once the QCP is located and the energy scale A de- 
termined as a function of ro — r 0ci we compute the gap 
m(T) and the pressure P(T), and deduce the universal 
scaling functions F^(x) and Fn{x) [Eqs. (4,3)]. To en- 
sure that we are in the universal (critical) regime, we 
solve the NPRG equations for various values of the ul- 
traviolet cutoff A, interaction strength uq or temperature 
T, and verify that the final results for Fn and Fn remain 
unchanged (Figs. 2 and 3). Only at sufficiently low tem- 
peratures and close enough to the QCP (T, |A| <C ck G ) 
do the universal scaling forms (3,4) hold. 

Figure 4 shows the universal scaling functions Fpj and 
Fn for various values of N (Table I shows the univer- 
sal ratio p s /(N\A\)). In the limit N — > 00, the trun- 
cated LPA' slightly differs from the exact result for the 
excitation gap m(T) but turns out to be extremely accu- 
rate for the computation of the pressure P(T) (the LPA' 
would be exact without the truncation of Uk(p) [36]). For 
smaller values of N, Fn and Fn differ significantly from 
the N — > 00 limit. While the large- AT result remains a 
good approximation in the quantum disordered regime, 
it becomes inaccurate in the quantum critical and renor- 
malized classical regimes. In particular, it misses the 
nonmonotonous behavior of F^{x) in the quantum crit- 
ical regime (|x| < 1) for N < 10. The possibility of such 
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x 

FIG. 5. Same as Fig. 4 but for negative values of x only. The 
dots show the large- N result rescaled by the factor (A — 1) /N . 



N 1000 10 8 6 4 3 2 

Cn/N to 0(1/N) 0.800 0.767 0.758 0.744 0.716 0.689 0.633 
Cn/N (NPRG) 0.812 0.796 0.793 0.788 0.781 0.775 0.767 



TABLE II. C N /N as obtained from the NPRG and the large- 
N approach [Eq. (47)]. 

a nonmonotonous behavior is discussed in Ref. [52]. 

In the renormalized classical regime, it is possible to 
reinterpret the large- N result so that it becomes consis- 
tent with the NPRG approach even for small values of N. 
Since the correlation length £ is exponentially large, we 
expect the thermodynamics to be dominated by the N— 1 
modes corresponding to transverse fluctuations to the lo- 
cal order. In the NPRG approach, these modes show up 
as Goldstone modes as long as po,k > (i.e. k > 
Since N — 1 is identified with N in the large- N approach, 
the latter overestimates the pressure, and therefore the 
scaling function Fn, by a factor N/ (N — l). In Fig. 5, we 
show that the large- N result, when rescaled by a factor 
(N — l)/N, is consistent with the NPRG approach. This 
shows that in the renormalized classical regime 

P(T)~(JV-1)^5' ( 45 ) 

which is nothing but the pressure of N — 1 free bosonic 
modes with dispersion uj — c|q| [53]. The very small 
excitation gap of the transverse fluctuations (m <C T) 
does not influence the thermodynamics. 

Of particular interest is the temperature variation of 
the pressure at the QCP (A = 0). Following Ref. [28], 
we express the pressure as 

P(T) = P(0) + C -^C N T ^ (46) 
27T cr 

for A = 0, where C N = N2ttF n (0)/((3). In the large-N 
limit [14], 

^ 4 _ 0,3344 /J_\ 
N 5 N \N 2 J { ' 



N = oo" 
1000 
10 



I 



K 




FIG. 6. Universal scaling function 3J-n{x) — xj-'(x) for the 
entropy per unit volume [Eq. (48)]. 



(the leading-order term 4/5 is given by Eq. (28)). The 
agreement between the 0(1/ N) result and the NPRG 
one rapidly deteriorates for N < 10 (Table II). 

The entropy per unit volume is equal to the tempera- 
ture derivative dP/dT of the pressure, 



S(T) 
V 



N- 



T 2 



T 



T 



T 



(48) 



Up to the factor NT 2 /c 2 , it is entirely determined by 
the universal scaling function 3Tn(x) — xT' N (x). The 
latter is nonmonotonous in the quantum critical regime 
(Fig. 6). 



C. RG flows 

In this section, we qualitatively discuss the RG flows in 
the various regimes of the phase diagram in the vicinity 
of the QCP for N > 3 (Fig. 1). We use the dimensionless 
variables 



Po, k = k L - d (Z A , k V A „ 
S k = (Z Ayk k 2 )~ l 8 k . 



a/2 



P0,k, 



(49) 



A fe = k^Z-fvxT 
and the corresponding RG equations 
dtPoM = (1 - d ~ W 2 - fj k /2)p~o,k 



3 f 



d t Sk 
d t ~\k 



— t, — h.t ii Po,k > 0, 



(%- 2)4 + y(^V + 2)4,i if Po.k 
(d-3 + 3 % /2 + %/2)A fc 
-Ag[9J fclU (0) + (iV-l)J fc> tt(0)], 



where 



m = 2\lp 0:k — [J kilt (p) + Jk,tl(y nljj=0 , 
Q 

fj k = 2X k p 0>k ^2 [Jk,it (p) + Jk,ti(P)]\p =0 , 



(50) 



(51) 
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FIG. 7. RG flows in the quantum critical regime (d — 2 and 
N = 3). The vertical dash-dotted line shows the thermal 
momentum scale kr = Ae -10 . [A = Co = 1, Uo = 27.6, 
r 0c ^ 0.065355.1 



with p = (p/k,Q n ), y = p 2 /fc 2 , Q n = uj n /c k k = 2nT k n 
(T k = T/c k k). The dimensionless threshold functions 
Ik, a and J k , a p(P) are defined in Appendix Bl. For the 
sake of generality, we consider an arbitrary space dimen- 
sion d. 

In the zero-temperature limit, using the results of Ap- 
pendix B 2 for the threshold functions, we recover the 
flow equations of the (d+l)-dimensional (classical) O(N) 
model in the LPA'. At the QCP (ro = ro c ), critical fluc- 
tuations develop below the Ginzburg momentum scale 
kQ. In the following, we discuss only the universal part 
of the flow k <C kQ. Deviations from criticality are 
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characterized by two momentum scales. The first one, 
&A = I A|/ c o, is associated to the detuning from the QCP. 
In the T = disordered phase, k^ 1 is nothing but the 
correlation length. In the T = ordered phase, k& ~ kj 
is related to the Josephson momentum scale kj = p s /cq. 
The latter separates the critical regime kj -C k <C &g 
from the Goldstone regime k <C fcj dominated by the 
Goldstone modes. The second characteristic momentum 
scale is the thermal scale kr = 2ttT/cq associated to the 
crossover between the quantum (k ^> kx) and classical 
(k -c hr) regimes. The three regimes of the phase dia- 
gram (Fig. 1) are defined by kr 3> &a (quantum critical), 
hr <C k/\ and r$ > rg c (quantum disordered), kx -C &a 
and r < ro c (renormalized classical). 



1. Quantum critical regime 

The RG flow in the quantum critical regime is shown 
in Fig. 7 for N — 3. The parameters of the microscopic 
action (1) are chosen such that the initial value A of the 
momentum cutoff is of the order of the Ginzburg scale 
ka- At the QCP (ro = ro c ) and for T = 0, we observe 
plateaus characteristic of critical behavior: pa j. ~ Pcrif 



A! 



and rjk = rjk = r\ (with rj the anomalous 



dimension at the three-dimensional Wilson-Fisher fixed 
point). At finite temperatures, the flow is modified when 
k becomes smaller than the thermal scale kx'- Po,k and 
rjkifjk rapidly vanish while Xk diverges; the (dimension- 
ful) order parameter po,fe vanishes and rat = y/SkJVAjc 
takes a finite value (indicating that the system is in a 
disordered phase) 



2. Quantum disordered regime 

Figure 8 shows the flow in the quantum disordered 
regime. At T = 0, the critical flow terminates at k ~ &a- 
For k < Po,k and rjk = fjk vanish while A^ di- 

verges; the (dimensionful) order parameter p ,fc vanishes 
and rrik = \f5k/VA,k takes a finite value. As expected, a 
finite temperature has hardly any effect on the flow when 
kx <C &a- Only for fcy ~ k& (i.e. near the crossover to 
the quantum critical regime) do we observe a modifica- 
tion of the T = flow. 



3. Renormalized classical regime 

The flow in the renormalized classical regime is shown 
in Fig. 9. The critical regime terminates at the Joseph- 
son scale kj. For k <§C kj, longitudinal fluctuations are 
suppressed (the (dimcnsionless) mass 2Xkf>o,k of the lon- 
gitudinal mode is much larger than unity) and the flow 
is dominated by the Goldstone modes. The anomalous 
dimensions and rjk then nearly vanish and the flow 
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FIG. 8. Same as Fig. 7 but in quantum disordered regime 
(ro = r-()c(l — 10~ 6 )). The dotted and dash-dotted vertical 
lines show the momentum scales &a an d fcy = Ae -12 , respec- 
tively. 



equations simplify into [54] 

Vk + Vk 



dtPo,k 



1 — d ■ 



P0,k 



N-l i 

2 KU (52) 



9 t A fc = (d - 3)A fe - X 2 k (N - 1) J fc ,tt(0), 

where the threshold functions are given in Appendix B 3. 
Po.k becomes very large in the renormalized classical 
regime. For k smaller than the inverse correlation length 
Po,k will ultimately vanish, but £ being exponen- 
tially large this is not seen in Fig. 9. In Sec. HID we 
discuss in more detail the behavior of p>o,fc and make the 
connection with the quantum NLcrM. 
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Ae" 



respec- 



For Ut *C k <C kj (quantum Goldstone regime), we 
can take the T = limit of the threshold functions (Ap- 
pendix B3a). We then find the fixed-point value 



A* 



d-3 



4?r 3/2 



(iV-l)J fc ,tt(0) (iV-l)(2- v / 2)' 



(53) 



where the last value is obtained with the exponential cut- 
off r(Y) = l/(e Y - 1) [Eq. (34)] and for d = 2. This 
fixed-point value shows up as a plateau A& ~ A* for 
kx <C k -C kj L which should not be confused with the 
plateau A& ~ A* rit corresponding to the critical regime 
kj <C k <C ka (Fig. 9). As discussed in detail in Ref. [25], 
the constant value Afc ~ A* and the diverging p ^ — k x ~ d 
(corresponding to a constant value of the order parame- 
ter po,fc) imply a vanishing of the longitudinal propaga- 
tor in the infrared limit: G\{p) ~ l/^ 2 , + c 2 p 2 )( 3_d '/ 2 
for k T <C |p|, |w n |/c <C fcj and d < 3 (the vanishing is 
logarithmic for d = 3). 

In the classical Goldstone regime k -C fop, the flow 
is dominated by classical (w n = 0) transverse fluctu- 
ations. As a result, the threshold function Jj, tt be- 
comes proportional to 71- (Appendix B 3 b) and A& van- 



ishes linearly with k. Since only the classical component 
4>{r) = 0(r, iw„ = 0) of the field matters, the effective 
action (37) becomes 



r c fe >] = pfd d 



J A,k 



At 



2 (V0r+y(p-^rj> (54) 



for po.fc > 0. Rescaling the field 
the usual form 



Tc(>, we obtain 



r£[0] = J d d r{ 



of the effective action for a classical model in the LPA', 
with the coupling constant A^. 1 = TX k - The appropriate 
dimensionless variable 



\cl _ \ cl 7-2 ld-i 
A k — A k Zj A,k K 

satisfies the RG equation 

d t ~Xf = (d - — (N — 1) J, c ; tt (0)(A c fe 1 ) 2 
with the threshold function 

Jk]tt — T k 1 JkM- 
This equation admits the fixed-point value 
d - 4 47T 



A 



(tf-l)J*(0) (JV-l)ln2' 



(56) 
(57) 
(58) 

(59) 



where the last value is obtained with the exponential cut- 
off and for d = 2. 



D. Goldstone regime and NLuM 

In the Goldstone regime, the behavior of the system is 
governed by the Goldstone modes and we expect a de- 
scription based on an effective NLctM to be possible. To 
identify the coupling constant of the effective NLctM [55] , 
we consider the following microscopic action 



%>] = / 
JO 



<h I d d r^{V V ) 2 
A 



V A 



+ ^(p-Po) 2 



(60) 



(p = cp 2 /2), which is analog to the ansatz (37,38) for the 
effective action r^[0]. This action can be written in the 
dimensionless form 



S{(p] = / 
Jo 



df 



•'■< \^ f Cp) 2 + \{dT$) 2 + ^{~p-~Pa,kf 



(61) 

where f = fcr, f = J Za/Va^^ and (p is defined in the 
same way as 4> in Eq. (49). Rescaling the field (p — > 
V / 2po<p, we then obtain 



S[<p] = Po df 
Ja 



d«A ( v^) 2 +(a^) 2 + A |°(^ 2 -i) 2 } 

(62) 
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In the limit Apo ^> 1, the last term in (62) imposes the 
constraint ip 2 = 1, and we obtain a quantum NLcrM with 
dimensionless coupling constant g = l/2p . 

In the NPRG approach, the coupling constant g k = 
l/2po,fe satisfies the RG equation 

dth = -(l-d- 2*±W^ ~ gk + ( N - i)/^, (63) 

which can be deduced from (52). Equation (63) should 
be considered together with the RG equation of the di- 
mensionless temperature T k = T/ckk, 

d t f k = - (l - J A-Hh^ f k . (64) 

Following Chakravarty et al. [18], we consider the cou- 
pling constants g k and t k = g k T k (rather than g k and 
fk), with 

d t t k = -(2 - d - m)h + (N — l)I k ,t9hh- (65) 

In the quantum Goldstone regime fc 3> fey, the system 
is effectively in the zero-temperature limit since C 1. 
Let us first consider the theta cutoff function [56] 

Rk(q) = z A , k (V - q 2 - <f\ e (V - q 2 - ^) ■ 

(66) 

In that case, one has 2rj k po,k = ~Ik,t (sec Ap- 
pendix B 3 c) , so that g k satisfies the flow equation 

d t9k = -{l-d)g k -2^\(N -2)~gl (67) 
a + 1 

We have used rj k — fj k (which holds for T k — !• 0) and 
J M = -2K d+1 /(d + 1) (with K d (27r) d = 4v d (2n) d the 
surface of the d-dimensional unit sphere). Equation (67) 
is nothing but the flow equation of the coupling constant 
in the d + 1-dimensional classical NLcrM [18, 57, 58]. It 
agrees with the result of Chakravarty et al. to order 
g\ [18]. In particular, we find that the 0(gj.) term van- 
ishes for the 0(2) model [N = 2) as expected. Note 
however that the coefficient of (N — 2)g 2 depends on the 
RG scheme and may therefore differ from the result of 
Ref. [18] obtained with a sharp cutoff. For an arbitrary 
function R k (q), the equality 2rj k po tk = —I k ± is in general 
violated, which does not allow us to recover the factor 
AT— 2 in the RG equation dtg k . For example, with the ex- 
ponential cutoff r(Y) = l/{e Y — 1), —2r]kPo,k/Ik,t — 1-11 
for d = 2. This failure is clearly an artifact of the LPA', 
which neglects the field dependence of Z A , k and Va,h- 
In the regime where Eq. (67) is valid, the C(ff 2 ) term is 
small and the LPA' remains nevertheless a very good ap- 
proximation to the flow equation of g k . In one dimension, 
d t g k is independent of the choice of the cutoff function 
R k to 0(g\), in agreement with the fact that the beta 
function of the classical two-dimensional NLtrM is uni- 
versal (i.e. independent of the RG scheme) to two-loop 
order [59]. 



There is a quantum classical crossover when k ~ kr, 
and for k <C fcr the system is governed by an ef- 
fective classical NLcrM with coupling constant t k [18]. 
With the theta cutoff (66), using 2r] k po. k = —I k ,t and 
/fc >t = —2{K d /d)i k /g k (Appendix B3c), we recover the 
beta function 

d t t k = -(2-d)t k -2^(N-2)t 2 k (68) 

of the (i-dimensional classical NLcrM [18]. As discussed 
above, for an arbitrary cutoff function R k (q), the coeffi- 
cient N—2is not exactly reproduced except for d = 2 due 
to the universality of the beta function d{t k to two-loop 
order in that case. 

E. 0(2) model: BKT transition temperature 

In the 0(2) model, there is a finite-temperature BKT 
transition for ro < tq c . For the classical 0(2) model, the 
NPRG reproduces most of the universal properties of the 
BKT transition [49, 50]. In particular one finds a value 
Pq of the dimensionless order parameter (the spin-wave 
"stiffness") such that the beta function /3(po k ) = d t po yk 
nearly vanishes for po ik > /5q. This implies the existence 
of a line of quasi-fixed points and enables to identify 
a low-temperature phase (T < Tbkt) where the run- 
ning of the stiffness Po, k , after a transient regime, be- 
comes very slow, implying a very large (although not 
strictly infinite as expected in the low-temperature phase 
of the BKT transition) correlation length £. In this low- 
temperature phase, the anomalous dimension r] k depends 
on the (slowly varying) stiffness po,fe- It takes its largest 
value ~ 1/4 when the RG flow crosses over to the disor- 
dered (long-distance) regime (for /5o,fc ~ Pg an d k ~ f" 1 ); 
and is then rapidly suppressed as po tk further decreases. 
On the other hand, the beta function is well approxi- 
mated by P(p ,k) = const x {p* - p ,k) 3/2 for p ,k < Pq, 
and the essential scaling £ - e con S t/(T-T BK T) 1/2 f tnc 
correlation length above the BKT transition tempera- 
ture ?bkt is reproduced [50]. Thus, although the NPRG 
approach does not yield a low-temperature phase with 
an infinite correlation length, it nevertheless allows us 
to estimate the BKT transition temperature from the 
value of /5q . A reasonable estimate of the BKT transition 
temperature in the two-dimensional XY model has been 
obtained using the NPRG [60]. The same method has 
been used to determine Tbkt in a two-dimensional Bose 
gas [51] in very good agreement with Monte Carlo simula- 
tions [61, 62]. We refer to Ref. [51] for more details about 
the determination of the BKT transition temperature in 
the NPRG approach. 

The BKT transition temperature corresponds to an es- 
sential singularity in the scaling function ^(x). Since 
is universal, the ratios Tbkt/|A| and Tbkt/Ps are also 
universal in the vicinity of the QCP (recall that p s = 
p s (T = 0) is the stiffness in the zero-temperature ordered 
phase). The NPRG approach predicts Tbkt/Ps — 1-59 
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FIG. 10. BKT transition temperature Tbkt vs the zero- 
temperature stiffness p s for various values of Mo [A = 100 
and co = 1]. 

with the exponential cutoff (Fig. 10) and Tbkt/ Ps — 1-5 
with a theta cutoff R k {q) = Z A , k {k 2 - q 2 )8(fc 2 - q 2 ) 
acting only on momenta. On the other hand, the ratio 
Tbkt / PsiT^rj,) = t/2 is universal anywhere on the tran- 
sition line, where Ps^bkt) denotes the stiffness jump at 
the transition [63]. While our determination of Tbkt 
is not precise enough to yield an accurate estimate of 
Tbkt/Ps, the latter is close to 7r/2, which implies that 
p s (Tg KT ) is only slightly reduced with respect to the 
zero-temperature stiffness p s . 

More generally, in the low-temperature phase near 
the QCP we can write the stiffness in the scaling from 
p s (T) = p s J{T/ Ps) with J[x) a universal scaling func- 
tion satisfying J'(O) = 1. The weak suppression of p s 
by thermal fluctuations for T < Tbkt implies that J{x) 
remains close to unity for x < 7rp s (Tg KT )/2/9 s . 



where the physics is dominated by the N — 1 Goldstone 
modes of zero-temperature broken-symmetry phase, the 
NPRG equations reproduce those of the quantum O(N) 
NLcrM [ ]. 

In the quantum 0(2) model, the ratio between 
the BKT transition temperature Tbkt and the zero- 
temperature stiffness p s (0) is universal near the QCP. 
The NPRG results show that Tbkt/Ps(0) is close to the 
universal ratio Tbkt / Ps(T^ KT ) = tt/2, implying that the 
stiffness Ps{Tbkt) a ^ the transition is only slightly re- 
duced with respect to p s (0). 

The superfluid-Mott-insulator (at constant density) of 
a Bose gas in an optical lattice provides us with a well 
controlled experimental realization of a relativistic QCP 
with a two-component (complex) field. Recent experi- 
ments have shown that it should be possible in the near 
future to observe quantum criticality associated with this 
QCP [64]. A measure of the temperature dependence 
of the pressure in the quantum critical regime would 
give an experimental estimate of the universal number 
^2(0) and enables a comparison with our theoretical re- 
sult J" 2 (0) ~ 0.147 (i.e. C 2 /2 ~ 0.767). Whether the full 
scaling function J-2(x) can be determined in the present 
experimental conditions requires a detailed study of the 
relativistic 0(2) QCP in the Bose-Hubbard model which 
will be reported elsewhere. 
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IV. CONCLUSION 

Using a NPRG approach, we have obtained the univer- 
sal function J-jv which determines the scaling form of the 
pressure near a relativistic QCP with O(N) symmetry 
[Eq. (3)]. For N < 10, the results are in strong disagree- 
ment with the large- N approach both in the renormalized 
classical and quantum critical regimes. If the large-A^ 
approach is properly interpreted, its results in the renor- 
malized classical regime can be reconciled with those of 
the NPRG approach. It fails however to describe the 
nonmonotonous behavior of the scaling function J-jy(x) 
in the quantum critical regime (\x\ < 1) as predicted by 
the NPRG approach. A similar nonmonotonous behavior 
is observed in the scaling function of the entropy. 

We have also shown how the NPRG allows us to ob- 
tain a complete picture of the quantum O(N) model in 
the vicinity of the zero-temperature QCP when N > 3. 
The characteristic momentum scales kx and &a, asso- 
ciated with temperature and detuning from the QCP, 
show up very clearly and yield distinctive RG flows in the 
quantum critical, quantum disordered and renormalized 
classical regimes. In the renormalized classical regime, 



Appendix A: Calculation of Trlngr 
In this appendix, we compute 

^Trlng- 1 = j m ( q 2 + Jjc? + m 2 /c 2 ). (Al) 

Since Trlng -1 is divergent, we subtract an infinite con- 
stant and consider 

Dim 1 ) = [ ln(q 2 + u?Jc 2 + m 2 /c 2 ) 

J q 

- [ [ Hq 2 +oj 2 /c 2 ). (A2) 

It is convenient to write D(m 2 ) — Do(m 2 ) + Di(m 2 ), 
where 

D (m 2 )= f f [ln(q 2 + W 2 /c 2 +m 2 / C 2 ) 

J q J to 

-ln(q 2 + W 2 /c 2 )], (A3) 
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and 



Di(m 2 ) = / ln(q 2 + oj 2 Jc 2 + m 2 /c 2 ) 

Jq 

ln(q 2 + lo 2 /c 2 + m 2 /c 2 ). (A4) 



Note that Di(m 2 ) vanishes at zero temperature. 



1. D (m 2 ) 



Using D o (0) = and 



D' (m 2 ) = 



1 



6r 0c m 



we obtain 



q , u uj 2 + c 2 q 2 + m 2 u c 2 4ttc 2 ' 

(A5) 



D (m 2 ) 



2 , 6ro c m 2 m 3 



UqC 2 6irc 2 ' 



(A6) 



Given that 



(a > 0), we deduce 



2. £>i(m 2 



da; 1 



1 1 



2?r w 2 + a 2 a C 30 - 1 



(A7) 



oxp ^ / 9y / c 2 q J - in- j 1 



Performing the momentum integral, we then obtain 
^A ? [- /3(c A-m)+ln(e^-l) 



(A8) 



D[(m 2 ) 



- ln(e^ m - 1)] 



(A9) 



for T < cA. From 
L>i(m 2 ) 



1 



ln(l 



2tt/3c 2 



2tt/3c 2 ^ fc 



(A10) 



we deduce 
L»i(m 2 ) = C* 

= C 



1 ^ //3m 1 \ 
tt/3 3 c 2 ^ V fc2 fc 3 / 

^[^Li 2 (e-^) + Li 3 (e-^)] 



(All) 



where lA s {z) is a polylogarithm [Eq. (27)]. The integra- 
tion constant C is fixed by requiring lim m 2_ > , 00 Di(m 2 ) = 
0. Since Li s (z) ~ z for |z| — » 0, this gives C = 0. 
3. £>(m 2 ) 

From Eqs. (A6,All), we finally obtain 



D(m 2 ) 



6r 0c m 2 



uqc 2 6irc 2 



7T/3 3 C 2 



[ / 3mLi 2 ( e ~' 3 ™) + Li 3 (e"' 3m )] . (A12) 



Using the same method, we can compute D{m 2 ) in the 
fermionic case (which amounts to replacing the bosonic 
Matsubara frequencies uj n by fermionic ones). We have 
verified that we then reproduce the result of Ref. [65] ob- 
tained by a different method. The fermionic result differs 
from the bosonic one only by the sign of the argument of 
the polylogarithm functions. 



Appendix B: Dimensionless threshold functions 

1. Definition 

The dimensionless threshold functions [Eq. (51)] are 
defined by 



du). 



p=0 



p=0 



2v d / y d/2 ~ 1 [r] k Yr + 2Y 2 r' + (fj k - m ){r + yr> 2 ]!" 2 , 



JkMO) = 2v d / y^-^mYr + 2Y 2 r' + (fj k - Vk )(r + Yr')u> 2 n ]A- 2 Af, 



= - 4 !f / y d/2 {2[ Vk Yr + 2Y 2 r' + (f ]k - Vk )(r + Yr')Q 2 n ]A' a A- 



- [%r + (% + A)Yr' + 2Y 2 r" + (% - m )(2r' + Yr")u 2 n ]A~ 2 }A' p Aj 2 , 
2v d f i// 2 - 1 [ Vk Yr + 2Y 2 r' + (fj k - %)(r + Yr')Cj 2 n ] {A\ - A 2 A )A- 2 Az\ 



(Bl) 
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where 



3. Goldstone regime 



A t =Y(l + r)+5 k , Ai = A t + 2\ k p^ k , 
A[ = A[ = 1 + r + Yr', 

and 

A 1 =2u> n (l + r + Yr'), 

A 2 = 1 + r + Yr' + 2w 2 (2/ + Yr"). 

We use the notations vj 1 = 2 d+1 ir d / 2 T(d/2), Y = y 
r = r(Y), r' = r'(Y), r" = r"(Y) and 



(B2) 



(B3) 



In the Goldstone regime 2X k po^ k 1, longitudinal 
fluctuations are subleading with respect to the transverse 
one. This yields I k i = Jfe n = 0, 



h,x = 4w d / y 
J§ 

Jk,tt = 4u<j / y 
Jp 



d/2-l Y 2 r , A -2^ 



d/2-ly2 7 ,/^-3 ) 



(BIO) 



and 



f^nj^rdy. (B4) JLj.^fi = V A f M 2r > + Yr") 1 + r + Yr ' 



2. Zero-temperature limit 

For T = 0, Lorentz invariance implies that ^A.fe = Va 
and r\ k =r) k . Using 



««,/ dy/ 72 " 1 



da) 
2tt 



u d+ i / dYY^y^fiY) (B5) 
/o 



for any function f(Y) = f(y + w 2 ), we obtain 

h,a = ^d+i / r (d+1)/2 (%r + 2Fr')i- 2 , 



Jk. a p(P) = 2%i / dYY( d+1 V 2 ( Vk r + 2Yr')A~ 2 A-„\ 
Jo 

(B6) 



Equation (B5) implies 
duo 



dyy 
Using 

we obtain 
d 



d/2 



f(Y) = V ^ dYY (d+1 ^ 2 f(Y). 
2tt V d +2 Jo 

(B7) 



Vd+3 d v d+1 



V d +2 d+1 v d 



(B8) 



dy 



[■h,\t{p) + Jk,ti(p)] 



p=0 



= -8 



Vd+l 

d+1 



dY y( d +!)/ 2 (i + r + Yr')A- 2 A t 



x [Y{ Vk r + 2Yr'){\ + r + Yr'^AV 1 + A^ 1 ) 

- Vkr - (Vk + 4)YV - 2Y 2 r"\ . (B9) 

Equations (B6) and (B9) yield the known RG equations 
of the (d+ l)-dimensional (classical) O(N) model in the 
LPA'. 



S-Jk,ti(p) 
dy 



kPo,k 

2 v d 

d .!,'' (l+r)» 



d/2 



1 + r + Yr' 



y 

x [2r' 2 -(l + r)r" 



(Bll) 



to leading order (we have set r\ k = r\ k = in the rhs of 
Eqs. (B10,B11)). 



a. Quantum Goldstone regime kr <C k 

For k T -C k (i.e. T <C c k k), we can take the zero- 
temperature limit in (BIO, Bll), which gives 



Jo (l + r) 2 

Jk,tt=4v d+1 dYY^/ 



(B12) 



(l + rf 



and 



■i-JfetfiP) 

dy 



2 v d 



+i 



5=0 W^+^o 

1 + r + Fr' , 



p=0 



(1 + 


r) 2 


2 




\ 2 n 2 


d+1 


1 + r - 


VYr' 


(1 + 


r) 3 



d yy(«J-l)/2 

(2r' + Yr"), 



[2/ 2 - (1 + r)r"] 



(B13) 

Equations (B12,B13) can also be deduced from (B6) and 
(B9) in the limit 2X k po.k ^ 1 and with S k = 0. 



b. Classical Goldstone regime k <C kr 

For k -C fey, the Matsubara sums in (BIO, Bll) are 
dominated by the zero-frequency term w„ = 0, which 
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gives 



h,t = 4vdT k / dyy 



d/2-1 







(1 + r) 2 



(B14) 



Jfc.tt = 4%T fc / dy y d/2 - 2 . 

lo {l + r) s 



and 



ay 



2 ^ 

-J-k 



dyy' 



d/2-i l + r + yr' 



x (2/ + yr") 



2 Vd* A 00 , rf/2 l + r + ^ 



p=0 



KpIm d Jo (1 + r) 3 



[2r' 2 - (1 + r)r"\ 



(B15) 



c. Theta cutoff 

In this section, we show that with the theta cutoff (66) 
the relation 1r\kPo.k = ~Ik,t is satisfied in the Goldstone 
regime (to leading order in l/A^po fc). Equation (66) im- 
plies 



r(Y) - ^— -0(1 -y) 



r'(Y) 



y 

1 e(i_y)_izZ«(i-y), 



y-2 



y 



(B16) 



At zero temperature, from Eq. (B12) we deduce 

o Vd + 1 



'fc.t 



1 



(B17) 



Since r(l + r + yr') = r'(l + r + yr') = 0, Eqs. (B13) 
simplify into 



p=0 
Vd+l 



■K-Jk,tl(P) 

dy 



1 



dyy(d+ l)/2 r ,, l+r + >V 

(1 + r)^ 



The product r"(l + r + Yr') is ill-defined with the theta 
cutoff because of the derivative 8'(1 — Y) of the Dirac 
function in r" . To circumvent this difficulty, we integrate 
by part, 



S-Jk,it{P) 
dy 



d 



p=0 
Vd+l 



Jk,t\{p) 



Vo,fc 



1 



dy 

dY < r 
I dY 



p=0 

9 y(d+i)/2 



1 + r 



+ 



r /2 Q y(rf+3)/2 

2" ay "(TTTp" 



i 



(B19) 



where we have used y(l + r) = 1 when < y < 1. 
From (B17) and (B19), we deduce 2r]kf>o,k — —Ik,t- 

In the classical Goldstone regime, we retain only the 
zero-frequency term Co n — in the Matsubara sums. The 
calculation of ifc :t and r\k is similar to the T = limit 
with the (d+ 1) -dimensional integrals over Y replaced by 
e?-dimensional integrals over y. Again we find ^rj^pQ^ = 
-/fc.t, with 



h.t = 4udT fc / dyy 1 
lo 



2-i r '(y) 

[l+r(y)} 2 



-8ffc^. 
d 



(B20) 
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